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Measurements of anisotropic flow Fourier coefficients («„) for inclusive charged particles and 
identified hadrons , K^, p, and p produced at midrapidity in Cu+Au collisions at ,/sffff = 200 
GeV are presented. The data were collected in 2012 by the PHENIX experiment at the Relativistic 
Heavy Ion Collider (RHIC). The particle azimuthal distributions with respect to different order 
symmetry planes f&n, for n = 1, 2, and 3 are studied as a function of transverse momentum px 
over a broad range of collision centralities. Mass ordering, as expected from hydrodynamic flow, 
is observed for all three harmonics. The charged-particle results are compared to hydrodynamical 
and transport model calculations. We also compare these Cu-fAu results with those in Cu-fCu 
and Au-|-Au collisions at the same and find that the V 2 and us, as a function of transverse 

momentum, follow a common scaling with 1/{snNpl^^). 

PACS numbers: 25.75.Dw 


I. INTRODUCTION 

Measurements of azimuthal anisotropies of particle 
emission in relativistic heavy ion collisions have proven 
to be an essential tool in probing the properties of the 
quark gluon plasma (QGP) produced in such collisions. 
These anisotropies can be quantified [T] by the coeffi¬ 
cients Vn in the Fourier expansion of the particle dis¬ 
tributions with respect to symmetry planes of the same- 
order Atn that are determined on an event-by-event basis: 
dN/d(j) oc 1 -I- X]n=i cos(n((/> — where n is the 

order of the harmonic, (/> is the azimuthal angle of parti¬ 
cles of a given type, and is the azimuthal angle of the 
n* **'-order symmetry plane. Measurements of the second 
harmonic, which indicates the strength of the “elliptic 
flow”, led to the conclusion that the QGP produced at 
RHIG behaves as a nearly inviscid fluid HE]. In the last 
decade, significant effort, both experimentally and the¬ 
oretically, has gone towards quantifying the specific vis¬ 
cosity r]/s (shear viscosity over entropy density) of the 
produced QGP, as well as its temperature dependence. 

Elliptic flow is thought to arise from the initial spa¬ 
tial anisotropy in the nuclear overlap zone, which has a 
lenticular shape in off-center nucleus-nucleus {A-\-A) col¬ 
lisions. This spatial anisotropy is then converted to a 
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momentum-space anisotropy through the pressure gradi¬ 
ents in the expanding fluid. Measurements of V 2 have 
been performed in symmetric A-\-A collision for a vari¬ 
ety of collision energies and particle species as a function 
of transverse momentum, rapidity, and system size H- 
HSj. Various scaling properties have been explored with 
the goal of understanding the onset of QGP formation 
with center-of-mass energy and how its properties may 
vary. The elliptic flow scaled by the corresponding initial 
spatial eccentricity (£ 2 ) was found to follow a universal 
trend when plotted against the produced particle den¬ 
sity in the transverse plane [a El [IS] over a broad range 
of center-of-mass energies. In a more recent study m, 
PHENIX showed that the transverse particle density is 
proportional to the third root of the number of partici¬ 
pant nucleons and that scaling with {e 2 N^i^^) re¬ 

moves the remaining system-size dependencies at various 
center-of-mass energies. 

The first-order Eourier coefficient vi, which is a mea¬ 
sure of the strength of the “directed flow”, has also been 
studied in symmetric A-\-A collisions over a broad range 
of energies in E HMH]- Most studies focus on mea¬ 
surements of pT-integrated values of vi as a function of 
rapidity or pseudorapidity, and the slope of dvi/dy at 
midrapidity, which may yield information on the location 
of a first-order phase transition in the phase diagram of 
nuclear matter [T5]. In symmetric A-\-A collisions, if the 
nuclei are considered to be smooth spheres, vi is an odd 
function with respect to (pseudo)rapidity and vanishes at 
midrapidity, which is consistent with the pT-integrated 
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measurements. 

Indeed, when the nuclei are taken as smooth spheres, 
all odd harmonics should vanish at midrapidity. How¬ 
ever, event-by-event fluctuations in the initial geometry 
can lead to nonzero odd harmonics at midrapidity m- 
Sizable values for these harmonics have been measured 
at both RHIC (ua) pTIE5] and the Large Hadron Col¬ 
lider {v 3 and U 5 ) [53H1S]- Evidence for a small rapidity- 
even component of vi at midrapidity has also been ob¬ 
served |18j . The combined experimental information 
from odd and even flow harmonics provides much more 
stringent constraints on the theoretical models 
and the extracted QGP properties than measurements 
of elliptic flow alone. 

Despite the wealth of experimental data and theoreti¬ 
cal studies, uncertainties in the energy density deposition 
in the initial state of the heavy ion collisions remain a lim¬ 
iting factor in deducing the specific viscosity of the QGP. 
Asymmetric collision systems, such as Cu-|-Au, provide 
opportunities to study the effect of the initial geometry 
on the collective flow, particularly because odd harmonics 
may be enhanced at midrapidity beyond the fluctuation 
effects. 

In this paper, we present measurements of vi, V 2 , and 
V 3 of charged particles and identified hadrons tt^, 
p, and p produced at midrapidity in Gu-|-Au collisions at 
y'Sjyjy = 200 GeV. In Sec.|^we present the experimental 
details of the measurements, and the sources of system¬ 
atic uncertainties. The results of the measurements are 
presented in Sec. m In Sec. |HI G| we compare the flow 
results obtained in different collision systems and explore 
their scaling behavior, and in Sec. HI D we present com¬ 
parisons to theoretical calculations. Sec. |IV| summarizes 
our findings. 


II. EXPERIMENTAL DETAILS 




18.5 m= 60 ft 


FIG. 1. The PHENIX detector configuration for RHIC Run- 
12 data taking period 


netic calorimeter (EMGal). The global system includes 
the beam-beam counters (BBGs), zero degree calorime¬ 
ters (ZDGs) and the shower maximum detectors (SMDs). 
Below, we give a brief description of each of these detec¬ 
tor sub-systems and their role in the present analysis. 


The PHENIX experiment is designed for the study of 
nuclear matter in extreme conditions using a wide vari¬ 
ety of experimental observables. The detector, optimized 
for the high-multiplicity environment of ultra-relativistic 
heavy ion collisions, comprises two central-arm spectrom¬ 
eters (East and West), two muon spectrometers (at for¬ 
ward and backward rapidity), and a set of detectors used 
to determine the global properties of the collisions. Eig- 
ure[^shows a schematic diagram of the PHENIX detector 
for the data recorded in 2012. The upper drawing shows 
a beam-axis view of the two central spectrometer arms, 
covering the pseudorapidity region |? 7 | < 0.35. The lower 
drawing shows a side view of the two forward-rapidity 
muon arms (North and South) and the global detectors. 
A detailed description of the complete set of detectors is 
given in Ref. [55] . 

The analysis presented here employs the global de¬ 
tectors, drift chamber (DG), three layers of multi-wire 
proportional chambers (PGI, PG2, and PG3), the time- 
of-flight detectors (TOEE, TOEW), and the electromag- 


A. Global Detectors 

The BBGs are located at ±144 cm from the nominal in¬ 
teraction point along the beam line, cover 27r in azimuth, 
and span the pseudorapidity range 3.0 < \r]\ < 3.9. Each 
BBG comprises 64 Cerenkov telescopes, arranged radi¬ 
ally around the beam line. The BBGs provide the main 
interaction trigger for the experiment and are also used 
in the determination of the collision vertex position along 
the beam axis (z-vertex) with cr^ = 0.6 cm resolution and 
the centrality of the collisions. The event centrality class 
in Cu±Au collisions is determined as a percentile of the 
total charge measured in the BBC from both sides. The 
BBGs also provide the start time for the time-of-flight 
measurement with a timing resolution around at = 40 ps 
in central Cu±Au collisions [55]. 

The ZDGs [34] are hadronic calorimeters located for¬ 
ward and backward of the PHENIX detector, along the 
beam line. Each ZDC is subdivided into three identi¬ 
cal modules of two interaction lengths. They cover a 
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pseudorapidity range of |? 7 | > 6.5 and measure the en¬ 
ergy of spectator neutrons with an energy resolution of 
a{E)/E = 85%/-\/i?-|-9.1%. The SMDs [SJ are scintilla¬ 
tor strip hodoscopes located between the first and second 
ZDC modules, a location corresponding approximately 
to the maximum of the hadronic shower. The horizon¬ 
tal coordinate is sampled by seven scintillator strips of 
15 mm width, while the vertical coordinate is sampled 
by eight strips of 20 mm width. The active area of each 
SMD is 105 mm x 110 mm (horizontal x vertical di¬ 
mension). Scintillation light is delivered to a multichan¬ 
nel Hamamatsu PMT R5900-M16 by wavelength shift¬ 
ing fibers [Mj. A typical position resolution for SMD is 
^ 0.1-0.3 cm. 


B. Tracking and Particle Identification Detectors 

The charged-particle momentum is reconstructed us¬ 
ing the tracking system. This system comprises the DC, 
located outside an axially-symmetric magnetic field at a 
radial distance between 2.0 m and 2.4 m, followed by 
PCl-3. The pattern recognition in the DC is based on 
a combinatorial Hough transform [35] in the track bend 
plane. A track model based on a field-integral look-up ta¬ 
ble determines the charged-particle momentum, the path 
length to the time-of-flight detector, and a projection of 
the track to the outer detectors. 

The tracks are matched to hits registered in the PC3 
and the EMCal, thus reducing the contribution of tracks 
originating from decays and y-conversions. 

The primary particle identification detectors used in 
this analysis are the time-of-flight detectors. The dif¬ 
ferent detectors in the east and west arms, use different 
technologies (scintillators and MRPCs respectively) and 
have different time resolutions [M1I3Z1- The total timing 
resolutions (including the start time measurement from 
the BBC) are 130 ps and 95 ps for east and west, respec¬ 
tively. Pion, kaon, and (anti)proton tracks are identified 
with over 97% purity for px < 2 GeV/c [Ml HI] in both 
systems. For px between 2-3 GeV/c, the purity of pions 
and protons is about 95% and that of kaons is around 
90%. 


C. Anisotropic Flow Measurement Technique 

The present measurements use the event-plane 
method |39j to quantify the azimuthal anisotropies of 
the particles produced in Cu-|-Au collisions. The vi, V 2 , 
and U 3 Fourier coefficients are determined as a function 
of centrality and px for inclusive charged particles and 
identified hadrons tt^, AT^, p, and p (with charge signs 
combined). 

In the event-plane method, a measured event- 
plane direction is determined for every event 

and for each order n. The harmonic coefficients 
Vni'^n} = (cosn((() — 4'°'^®)) /Resjdin} are then mea¬ 

sured with respect to the event plane for each har¬ 
monic, where (j) is the azimuthal angle of the hadron and 
Resj'I'n} is the event-plane resolution. 

The collision geometry of a Cu-|-Au collision is shown 
in Fig. ^a) projected onto the reaction plane, and in 
Fig. I^b) projected onto the plane perpendicular to the 
beam axis. Figure [^a) shows direction of the projectile 
(Gu) and target (Au) spectators, which are bent away 
from the participant zone. There is an alternative pic¬ 
ture, in which the spectators are attracted towards the 
center of the system, as discussed in m- In this paper, 
we assume the former picture for the determination of the 
direction of the event-plane angle from the spectators. 

As shown in Fig. [^a), the Cu spectators fly along the 
the positive-rapidity direction (North) and the Au spec¬ 
tators go towards the negative-rapidity direction (South). 
The central position of the Au spectators is measured 
by the South SMD (SMDS) to determine the spectator 
plane gj charged and identified hadrons 

is measured with respect to 4 /Smds^ g^g indicated in Eq. 
Measurement with respect to the spectator plane is pre¬ 
ferred over the first order event-plane determined by the 
distribution of the produced particles, because the distri¬ 
bution of the spectators is less distorted by momentum 
conservation effects. 


= -(C0s(</.t,ack - 4 /SMDS)^/j^gg(^SMDS) 


Res(4'f“°®) 


(cos I 


(•^SMDS _ 


(cos (^SMDS _ ^SMDN^^^ggg ^^SMDS _ ^BBCSN^^ 

(C0S(^SMDN_^BBCSN^^ 


( 2 ) 


There is a negative sign in Eq. to keep the conven¬ 
tion in which the direction of projectile (Cu) spectators 
is positive. In Eq. ^ the resolution of 4 /Smds jg 
lated in 10% centrality intervals with the three subevent 
method 133 ED by combining the other Cu spectator 
plane from the North SMD (SMDN) and the first or¬ 


der participant event-plane measured by the combined 
South and North BBCs (BBCSN). However, this method 
for determining the resolution assumes a nonfluctuating 
nuclear-matter distribution. Event-by-event fluctuations 
in the initial energy density of the collision will cause 
the vi signal to be different with respect to 4 /Smds 
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FIG. 2. Sketch of a noncentral heavy-ion collision. See text for description of the figure. 
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FIG. 3. Panel (a) event-plane resolution as a function of centrality for the SMDS detectors. Panel (b) and (c) second- and 
third-order event-plane resolution. The BBC event-plane resolution is obtained from two sub-events and BBCS, BBCN, CNT 
from three sub-events as a function of centrality. 


T'fMDN .j-Q rapidity-symmetric component in the 
direct flow [32]. To cover this uncertainty, the resolution 
of is also calculated using the participant plane 

from either BBCS or BBCN and the differences are as¬ 
signed as a systematic uncertainty. 

The second ('I> 2 ) and third (T'a) order event planes are 
measured by the combination of BBCS and BBCN. To 
determine the second and third order event-plane reso¬ 
lution from the BBC, we first measure the second and 
third order event planes with the BBCS (Au-going side), 
BBCN (Cu-going side) and central arm tracks (CNT). 
The central-arm tracks are restricted to low pT (0.2 
< Pt < 2.0 GeV/c) to minimize the contribution from jet 
fragments. The second and third order event-plane res¬ 
olution of BBCS, BBCN, and CNT are calculated using 
three subevent methods with a combination of BBCS- 
BBCN-CNT. Then the second and third order event- 
plane resolutions of the BBC (including both BBCS and 
BBCN) are calculated with two subevent methods with 


a combination of BBC-CNT. 

The event-plane resolutions for different subsystems 
are shown in Fig.j^as a function of centrality. Panel (a) of 
Fig.ll shows the resolution of the first-order event plane 
as measured by the SMDS using three different meth¬ 
ods. The first method uses a three subevent combination 
SMDS-BBCSN-SMDN, shown with circles, the second 
method shown with open squares uses a three subevent 
combination SMDS-BBCS-SMDN, and the third method 
shown with open triangles uses the combination SMDS- 
BBCN-SMDN. The resolution of the second and third 
order event planes for BBC, BBCS, BBCN, and CNT 
are shown in panel (b) and (c) of Fig. respectively. 


D. Number of Participants and Eccentricity 

A Monte-Carlo Glauber simulation was used to esti¬ 
mate the average number of participating nucleons Apart 
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and the eccentricity 

\/ (r^ cos(n^))2 + (r^ sin(n0))2 

This simulation employed a Glauber model with a 
Woods-Saxon density profile and includes modeling of 
the BBC response [33JI33]. The eccentricity defined in 
Eq. is also known as the participant eccentricity ffpart 
and includes the effect of fluctuations from the initial 
participant geometry. Table summarizes fVpart and e„. 

E. Data set 

The measurements presented here use data from 
Cu+Au collisions at = 200 GeV collected by the 

PHENIX experiment at RHIC in 2012. Minimum bias 
events triggered by the BBC recorded within ± 30 cm 
from the nominal interaction point along the z-axis were 
used. The events were examined to ensure that stable 
performance is seen in the detectors used in the analy¬ 
sis, namely DC, PCS, TOE, BBC, and SMD. A total of 
3.6 X 10® events were analyzed. 

F. Systematic Uncertainties 

Tables [llHV] summarize the systematic uncertainties 
for the measurements of Vi, V 2 , and for inclusive and 
identified charged hadrons, which are categorized by the 
types: 

A point-to-point uncertainties uncorrelated between 
Pt bins; 

B pT-correlated, all points move in a correlated man¬ 
ner, but not by the same factor; 

C an overall normalization error in which all points 
move by the same multiplicative factor independent 
of Pt- 

Contributions to the uncertainties are from the follow¬ 
ing sources: 

1 . event-plane resolution correction, 

2 . event plane as measured using different detectors, 

3. Vn from background tracks, 

4. acceptance dependencies 

5. PID purity. 

The uncertainties from measurements of the event 
planes using different detectors are found to only weakly 
depend on px- For the measurement of Vi, the uncer¬ 
tainties are obtained by comparing the vi as measured 


with SMDS with alternately BBCN or BBCS used for 
resolution. For V 2 and fa, the uncertainties are obtained 
by comparing the V 2 and as measured by the BBCN 
and BBCS. For the vi measurement, for the 10%-20% 
centrality class we find a 20 % systematic uncertainty in¬ 
dependent of Pt- For the 40%-50% centrality class, we 
find a 12% systematic uncertainty. For V 2 , the system¬ 
atic uncertainty is less than 3% for the 0%-10% centrality 
range and increases to 4% for the centrality range 50%- 
60%. For f 3 , a 3% systematic uncertainty is found for 
0%-10% centrality, increasing to 7% for the 20%-30% 
centrality range. 

Background tracks that are not removed by the track¬ 
ing selections as described in Sec.H may influence the 
measured Vn- They can arise from particle decays, 7 - 
conversions, or false track reconstruction. We estimate 
the tracking background contribution by varying the 
width of the track-matching window in PC3 and com¬ 
paring the results with and without the EMCal matching 
cut. We find that the absolute uncertainty for vi is less 
than 5 x 10“^. For V 2 and V 3 , the change is less than 2%. 

Systematic uncertainties of acceptance were evaluated 
using different subsets of the detector such as DC and 
TOF in the east and west arms. Differences in the f„ 
measured using different arms may be caused by differ¬ 
ent detector alignment and performance. Maximum dif¬ 
ferences of order 3% and 10% were found for V 2 and V 3 
respectively. These uncertainties have centrality depen¬ 
dence and minimal px dependence. For vi, maximum 
absolute uncertainty of 3 x 10“^ is found. These uncer¬ 
tainties are detailed further in Tables El and mil 

An additional systematic uncertainty in v-n resulting 
from hadron misidentification is based on the PID purity 
estimates from the TOF detectors as discussed in Sec.II. 
Pion, kaon, and proton species purity is greater than 90% 
and the differences between their corresponding v-n is less 
than a factor of two. For V 2 and 113 , an additional uncer¬ 
tainty of 3% (type A) attributable to contamination from 
other species is found for particles with px < 2 GeV/c, 
5% for higher px pions and protons, and 10% for higher 
Px kaons. In the measurements of z;i, a common abso¬ 
lute uncertainty of 1 x 10 “^ is found for the three particle 
species for px < 2 GeV/c, and at higher px the uncer¬ 
tainties are 2x10“® for pions and 3x10“® for kaons and 
protons, respectively. The uncertainties due to particle 
identification are to be added in quadrature to the values 

listed in Tables El and ED 


III. RESULTS AND DISCUSSION 

A. Harmonic flow results from Cu-f Au collisions 

Figures 14^^ show the vi, V 2 , and V 3 results for charge- 
combined hadrons measured as a function of px in 
Cu-I-Au collisions at y/s^ = 200 GeV. Different cen¬ 
trality intervals are studied. The filled circles show the 
Vnipx) data, and the systematic uncertainties are shown 






TABLE I. Number of participants and the participant eccentricity (£ 2 , £ 3 ) from Monte-Carlo Glauber calculations for Au-|-Au, 
Cu-bCu and Cu+Au collisions at 200 GeV 


centrality 

bin 

-^part 

Au±Au 200 GeV 

£2 

£3 

Cu±Gu 200 GeV 

-^part ^2 

-^part 

Cu±Au 200 GeV 

£2 

£3 

0%-10% 

325.2 

0.103 

0.087 

98.2 

0.163 

177.2 

0.138 

0.130 


±3.3 

±0.003 

± 0.002 

±2.4 

±0.003 

±5.2 

±0.011 

±0.004 

10%-20% 

234.6 

0.200 

0.122 

73.6 

0.241 

132.4 

0.204 

0.161 


±4.7 

±0.005 

± 0.004 

±2.5 

±0.007 

±3.7 

±0.008 

±0.005 

20%-30% 

166.6 

0.284 

0.156 

53.0 

0.317 

95.1 

0.280 

0.208 


±5.4 

±0.006 

± 0.005 

±1.9 

±0.006 

±3.2 

±0.008 

±0.007 

30%-40% 

114.2 

0.356 

0.198 

37.3 

0.401 

65.7 

0.357 

0.266 


±4.4 

±0.006 

± 0.008 

±1.6 

±0.008 

±3.4 

±0.010 

±0.010 

40%-50% 

74.4 

0.422 

0.253 

25.4 

0.484 

43.3 

0.436 

0.332 


±3.8 

±0.006 

± 0.011 

±1.3 

±0.008 

±3.0 

±0.013 

±0.013 

50%-60% 

45.5 

0.491 

0.325 

16.7 

0.579 

26.8 

0.523 

0.412 


±3.3 

±0.005 

± 0.018 

±0.9 

±0.008 

±2.6 

±0.019 

±0.019 


0.02 

0 

^- 0.02 
-0.04 

012345012345012345012345 
(GeV/c) (GeV/c) p^ (GeV/c) p^ (GeV/c) 

FIG. 4. v\{jpt) for charged hadrons measured with respect to the Cu spectator neutrons at midrapidity in Cu-I-Au collisions 
yJ^NN ~ 200 GeV. Error bars show the statistical uncertainties, and shaded boxes indicate the systematic uncertainties. 
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TABLE 11. Systematic uncertainties in the ui measurements. 


Wi Uncertainty Sources 

10%-20% 

40%-50% 

Type 

Ui Event-plane 

20% 

12% 

C 

Background(absolute value) 

5xl0~^ 

bxlO""^ 

A 

Acceptance (absolute value) 

3xl0~^ 

2x10"^ 

C 


with the shaded boxes. 

The vi{pt) measurements shown in Fig. are per¬ 
formed with respect to the event plane determined by 
spectator neutrons from the Au nucleus. To align with 
previous conventions, we flip the sign so that it is effec¬ 
tively with respect to the spectator neutrons from the Cu 
nucleus, as noted in Sec. |II C[ In all centrality intervals, 
high pt particles at midrapidity move in the direction 
opposite of the Cu nucleus spectator neutrons, as indi- 


TABLE III. Systematic uncertainties given in percent on the 
V 2 and U 3 measurements. 


Un(n=2,3) Uncertainty Sources 

0%-10% 

to 

0 

1 

CO 

0 

Type 

V 2 Event-plane 

3% 

4% 

B 

Background 

2% 

2% 

A 

Acceptance 

2% 

3% 

C 

113 Event-plane 

3% 

7% 

B 

Background 

2% 

2% 

A 

Acceptance 

8% 

10% 

C 


cated by the negative vi values. Low px particles might 
then be expected to move in the opposite direction by 
conservation of momentum, and there is a hint of this 
effect though not beyond current systematic uncertain- 
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FIG. 5. V 2 {pt) for charged hadrons measured at midrapidity in Cu+Au collisions at = 200 GeV. Uncertainties are as 

in Fig.|4] 
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FIG. 6. V 3 {pt) for charged hadrons measured at midrapidity in Cu+Au collisions at = 200 GeV. Uncertainties are as 

in Fig. |4] 
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TABLE IV. Systematic uncertainties in the measured Vi for TABLE V. Systematic uncertainties in percent on the mea- 
identified particles. sured V 2 and V 3 for identified particles. 


species 

Pt < 2 GeV/c pt 

> 2 GeV/c Type 

species 

Pt <2 GeV/c 

Pt > 2 GeV/c 

Type 

pion (absolute value) 

1 x10"® 

2x10"® A 

pion 

3% 

5% 

A 

kaon (absolute value) 

1 x10"® 

3x10"® A 

kaon 

3% 

10 % 

A 

proton (absolute value) 

1 x10"® 

3x10"® A 

proton 

3% 

5% 

A 


ties. The Vi component is consistent with zero for < 1 maximum of the absolute value decreases from central to 
GeV/c and its absolute value increases at higher pt. The peripheral collisions. This is contrary to the centrality 
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dependence of V 2 where the values increase from the most 
central 0%-10% collisions, up to the 30%-40% centrality 
class. This trend in V 2 is expected from the initial ge¬ 
ometry, because the ellipticity of the participant zone £2 
(see Table |T| increases in the peripheral collisions. The 
V 2 {pt) values in the 30%-40%, 40%-50%, and 50%-60% 
Cu-l-Au centrality classes, shown in Fig. are consistent 
with each other, showing very little, if any centrality de¬ 
pendence. The V 2 and values are positive, as previ¬ 
ously observed in symmetric collisions systems. For all 
three harmonics, the magnitude of the signal increases 
with pt up to about px = 3 GeV/c, and then tends to 
decrease. This may indicate a change in the dominant 
production mechanism, e.g., an increasing contribution 
from jet fragments, or it may be due to the fact that 
higher px particles escape the fireball with fewer interac¬ 
tions. 

The V 3 component (Fig. has weak centrality depen¬ 
dence, a behavior which is similar in symmetric A+A col¬ 
lisions [211 [ 22 ] , where the triangular flow at midrapidity 
is completely driven by the event-by-event fluctuations of 
the interaction zone. These fluctuations are also present 
in the asymmetric Cu-|-Au collisions and are expected to 
play a similar role. In Sec. |III C| we compare the flow re¬ 
sults obtained in different collisions systems and explore 
their scaling behavior. 


B. Identified particle flow results 


Figures!^ andshow the particle-species dependence 
of V 2 and V 3 in Cu-|-Au collisions. Results are presented 
for charge-combined tt^, p, and p. The measured 
Vnipr) values are shown with points, and the shaded 
boxes represent the species-dependent type A systematic 
uncertainties. The type B and C systematic uncertain¬ 


ties shown in Table III are largely common for all particle 
species. For the odd harmonics, to improve the statisti¬ 
cal significance of the results the measurements for identi¬ 
fied particles are performed in a single centrality interval, 
namely 0%-30% for v^ijpx) and 10%-50% for vi{px)- 
There are two trends common to both n = 2,3 re¬ 
sults shown in Figs. and First, in the low-p^ region 
the anisotropy appears largest for the lightest hadron 
and smallest for the heaviest hadron. A similar mass 
ordering is also predicted by hydrodynamics, in which 
all particles are moving in a common velocity field. Sec¬ 
ond, for Pt > 2 GeV/c this mass dependence is reversed, 
such that the anisotropy is larger for the baryons than 
it is for mesons at the same px- These patterns have 
been observed previously in Vn measurements for identi¬ 
fied particles in Au-|-Au collisions at RHIC. The vi{px) 
values, presented in Fig. also show mass ordering, al¬ 
though these measurements have larger overall system¬ 
atic and statistical uncertainties than V 2 {px) and v^(px)- 
As in t he c ase of vi{px) for charged particles described 
in Sec. |ni] (Fig. i. we note that although the values of 
vi{px) for each species appear to be positive at low px, if 


the full systematic uncertainty of type B and C is taken 
into account, a definitive conclusion can not be drawn 
about the overall sign of the bulk directed flow. The 
mass dependence in the collective flow at the low-pT is 
a generic feature of hydrodynamical models. The depen¬ 
dence on valence quark number in the intermediate-pT 
region has been associated with the development of flow 
in the partonic phase of the fireball evolution and subse¬ 
quent hadronization by parton coalescence [45] . 


C. System size dependence 


It is interesting to compare the charged-hadron Vnipx) 
results for different collision systems measured in the 
same experiment at the same center-of-mass energy. 
PHENIX has previously studied anisotropic flow har¬ 
monics in symmetric Au-|-Au and Gu-|-Cu collisions at 
,/s^ = 200 GeV [m 121] . By varying the system size 
and the centrality selection, one can study the effects of 
the initial geometry on the observed flow coefficients. We 
will first compare the results obtained in different colli¬ 
sion systems for the same centrality selections, and then 
explore possible scaling behaviors. 

In Fig. 10 the V 2 {px) coefficients are compared for six 
different centrality selections. We observe that in each 
centrality class at a given px the values measured in 
Cu-|-Au collisions are always between those measured in 
Cu-|-Gu and Au-|-Au collisions. In all centrality classes 
chosen, the Cu-|-Gu system has larger elliptic eccentricity 
than both Cu-|-Au and Au-|-Au collisions. However, ex¬ 
cept in the most central 0 %- 10 % collisions, the measured 
V 2 {pt) values are not ordered according to the magnitude 
of £2 in the different systems listed in Table |lj To further 
investigate this, in Fig. 11 we scale the V 2 {px) values 
in each collision system with their respective participant 
eccentricity £ 2 - The resulting V 2 (px )/£2 are ordered by 
system size, but this scaling does not lead to a universal 
behavior. 

In Ref. [12], PHENIX compared measurements in 
Cu-|-Gu and Au-l-Au collisions for different center-of- 
mass energies and centrality selections and found that the 

1/3 

V 2 values obey common empirical scaling with e 2 A^part- 

1 /3 

The motivation for introducing the factor is that 

under the assumption that Vpart is proportional to the 
volume of the fireball, is a quantity proportional to 

a length scale, and therefore may account for the system- 
size dependence of the V 2 values. In Fig. we add to 
this comparison the results from the asymmetric Cu-|-Au 
collisions. This scaling brings the V 2 {px) results from 
the three collisions systems together across all centrality 
classes in this study. 

In Fig. 13 the vz{px) values are compared in Cu-|-Au 
and Au-|-Au collisions for events of the same centrality. 
Unlike in the V 2 {px) measurements, here the values of 
v^{j)t) are ordered according to the initial triangularities 
£3 listed in Table |Tj with the Gu-|-Au results being larger 









11 



p^(GeV/c) PT{GeV/c) P-r(GeV/c) 


FIG. 7. The second-order Fourier coefficients V 2 {pt) for charge-combined identified hadrons p, and p measured at 

midrapidity in Cu-|-Au collisions at ,^8^^ ~ 200 GeV for the centrality classes marked in each panel. The symbols represent 
the measured V 2 {pt) values, the error bars show the statistical uncertainties, and the shaded boxes indicate the systematic 
uncertainties from PID. The full systematic uncertainties, that are mostly common to all particle species are shown in Table |In| 



FIG. 8. The third-order Fourier coefficients v^ipr) for 
charge-combined identified hadrons tt^, p, and p mea¬ 
sured at midrapidity in Gu+Au collisions at = 200 

GeV for 0%-30% centrality. The symbols represent the mea¬ 
sured vz(pt) values, the error bars show the statistical uncer¬ 
tainties, and the shaded boxes indicate the systematic uncer¬ 
tainties from PID. The full systematic uncertainties, that are 
mostly common to all particle species are shown in Table [Till 
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FIG. 9. The first-order Fourier coefficients vi(pt) for charge- 
combined identified hadrons tt^, K^, p, and p measured at 
midrapidity in Cu-I-Au collisions at = 200 GeV for 

10%-50% centrality. The symbols represent the measured 
v\{pt) values with respect to the Gu spectator neutrons, the 
error bars show the statistical uncertainties, and the shaded 
boxes indicate the systematic uncertainties from PID. The 
full systematic uncertainties, that are mostly common to all 
particle species are shown in Table [III 
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(GeV/c) (GeV/c) p^ (GeV/c) 

FIG. 10. The second-order Fourier coefficients V 2 {pt) for charged hadrons measured at midrapidity in Cu-|-Au, Au-|-Au |12| . 
and Cu-I-Cu m collisions at = 200 GeV. In each panel, the V 2 {pt) coefficients are compared for the same centrality class, 

as marked in the figure. The symbols represent the measnred V 2 {pt) values, the error bars show the statistical uncertainties, 
and the shaded boxes indicate the systematic uncertainties. 



FIG. 11. Scaled second-order Fourier coefficients U 2 (pt )/£2 for charged hadrons measured at midrapidity in 

Cu-|-Au,Au-|-Au m, and Cu-|-Gu m collisions at = 200 GeV. In each panel, the V 2 {pt) values measured in the 

centrality classes marked in the hgure, are scaled by the average second-order participant eccentricity £2 in the initial state of 
the collisions as determined by a MG Glauber calculation described in the text. The symbols represent the scaled U 2 (pt )/£2 
values, the error bars show the statistical uncertainties, and the shaded boxes indicate the systematic uncertainties. 
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(GeV/c) (GeV/c) p^ (GeV/c) 

FIG. 12. Scaled second-order Fourier coefficients V2{pT)/{£2N^iyt) charged hadrons measured at midrapidity in 
Cu-l-Au,Au-|-Au m, and Cu-|-Cu m collisions at = 200 GeV. In each panel, the V 2 {pt) values measured in the 

centrality classes marked in the figure, are scaled by the average second-order participant eccentricity £2 in the initial state 
of the collisions as determined by a MG Glauber calculation described in the text, and the corresponding number of nucleon 
participants The symbols represent the scaled V2{pt)/{£ 2N^i^t) values, the error bars show the statistical uncertainties, 

and the shaded boxes indicate the systematic uncertainties. 



FIG. 13. The third-order Fourier coefficients V 2 ,{pt) for charged hadrons measured at midrapidity in Gu+Au and Au-fAu |22| 
collisions at = 200 GeV. In each panel, the v^ipr) coefficients are compared for the same centrality class, as marked in 

the figure. The symbols represent the measured W 3 (pt) values, the error bars show the statistical uncertainties, and the shaded 
boxes indicate the systematic uncertainties. 
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than the Au+Au ones. In particular, in the most central 
0%-10% collisions £3 in Cu+Au is about 50% larger than 
in Au+Au collisions, and a similar difference is observed 
in the vz{pt) values. In Fig. 14 the vz(pt) values are 
scaled by the initial £3 eccentricity. A good agreement 
between the different systems is observed at low px (< 2 
GeV/c), which indicates that the participant eccentrici¬ 
ties obtained in the Glauber model provide an adequate 
description of the fluctuating initial geometry. Addition- 
ally, we perform scaling with as was done for 

the V 2 {pt) measurements. The results of this scaling 
are shown in Fig. 15 In this case, the measurement in 
Cu+Au and Au+Au collisions are in better agreement 
at high pt, however at low px the v^ijpx )/values 
are systematically higher for the Cu+Au system. 


D. Theory comparisons 

1. Hydrodynamic calculations 


Predictions from 3D+1 viscous hydrodynamic calcula¬ 
tions are available [15]. At lowp^ (< 1-0 GeV/c) directed 
flow is predicted to be in the hemisphere of the Cu side, 
while for high px ( > 1.5 GeV/c) directed flow is pre¬ 
dicted to be in the hemisphere on the Au side. Further, 
the bulk directed flow component from integration over 
Px is predicted to be in the Cu-nucleus hemisphere. Due 
to the large systematic uncertainties and small value of 
Vi at small px^ we can not reliably determine the sign 
of the vi component at low px, or the sign of the bulk 
directed flow. At high px the measurement is in agree¬ 
ment with the directed flow being in the Au hemisphere, 
under the assumption that the spectator neutrons are de¬ 
flected outward from the interaction region and aligned 
with the impact parameter vector. Ref. [46] shows the 
vi with respect to the reaction plane (i.e. the impact pa¬ 
rameter vector) for 20%-30% central Cu+Au collisions 
including particles within \ri\ < 1 . 0 , and thus we cannot 
compare directly with our narrower rapidity selection. It 
is notable however, that the hydrodynamic results at px 
= 2 GeV/c reach vi « 5%, while the experimental data 
within |? 7 | < 0.35 are less than 2%. 

The predictions for elliptic and triangular flow as a 
function of px are compared to the data in Fig. 16 and 
Fig. Calculations with two different values of the 
specific viscosity rj/s = 0.08 and rjjs = 0.16 are shown. 
Our measurements in the 20%-30% centrality range are 
consistent with each of these values; for the most central 
0%-5% events, a value of ry/s = 0.08 is closer to the data. 


2. AMPT 

The A-Multiphase-Transport Model (AMPT) genera¬ 
tor gaiiHi has been established as a useful tool in the 
study of flow observables in heavy-ion collisions [45] . 
Therefore, it is of interest to compare the measured ui, U 2 , 


and V 3 as a function of px with the corresponding quan¬ 
tities calculated using the AMPT model. To that end, 
we used AMPT v2.21 with string melting turned on to 
generate approximately 2 million minimum bias Cu+Au 
events at = 200 GeV, setting the partonic cross 

section alternately to Cpart = 1-5 mb and 3.0 mb. In 
the default version of the model, initial conditions are 
generated using Monte Carlo Glauber with a gray disk 
approach to nucleon-nucleon interactions. However, in 
this study we utilize a modified black disk Glauber model 
with a fixed nucleon-nucleon inelastic cross section of 42 
mb, as used in Ref. [45] . 

Following the method of [25] , Fourier coefficients vi,V 2 , 
and V 3 are calculated for unidentified charged hadrons 
within |ry| < 0.35, with respect to the corresponding par¬ 
ticipant planes and 4 > 3 . These plane angles are 

computed for each event from the initial coordinates of 
nucleon participants with a Gaussian smearing of width 
a = 0.4 fm. 


The V 2 (px) and vs^px) results shown in Fig. 18 and 
Fig. are well reproduced by the model for px < 1 
GeV/c. The comparison with the data indicates that 
the 3.0 mb partonic cross section gives a better descrip¬ 
tion of the system dynamics. However, the calculation 
of vi and its comparison with experimental data is less 
straightforward. Because the experimentally measured 
Cu spectator neutron orientation is unknown, we calcu¬ 
late the vi values with respect to the impact parame¬ 


ter vector b pointing in the direction of the Cu nucleus 
as well as with respect to dti, the overlap region calcu¬ 
lated as previously described. Because the calculation is 
done in the participant center-of-mass frame, weighting 
all participants equally yields exactly £1 = 0 and hence 
no direction for dti. There are various suggestions in the 
literature for weighting with and r3 [Snilil], and in 
this study we choose to use r^. 

In addition, we have considered two different Monte 
Carlo Glauber initial conditions, one with black disk 
(BD) nucleons and one with gray disk (GD) nucleons, 
thus varying the diffuseness of the nucleon-nucleon inter¬ 
action radius. Figure [^ shows results for Cu+Au col¬ 
lisions within the 30%-40% centrality selection on the 


relative distribution of 4 '1 to 6 pointing the direction of 
the Cu nucleus. Panel (a) is for the BD case and Panel 
(b) the GD case. This small difference in the treatment 
of initial geometry completely re-orients the dtj vector. 
The lower panels show the AMPT midrapidity particle 
vi as a function of px relative to 4'i and b again the BD 
and GD implementation. It is interesting to note that in 
the GD case where the two results agree, the prediction 
is for low Px particles moving in the direction of the Au 
nucleus and the high px particles in the direction of the 
Cu nucleus (opposite to the previously discussed hydro- 
dynamic prediction). 
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FIG. 14. Scaled third-order Fourier coefficients V 3 {pt)/S 3 for charged hadrons measured at midrapidity in Cu-|-Au and 
Au-I-Au m collisions at = 200 GeV. In each panel, the vzipr) values measured in the centrality classes marked in the 

figure, are scaled by the average third-order participant eccentricity £3 in the initial state of the collisions as determined by a 
MG Glauber calculation described in the text. The symbols represent the scaled V 3 {pt)/s 3 values, and the error bars show the 
statistical uncertainties. The shaded boxes indicate the systematic uncertainties in the Cu-|-Au measurements, and the lines 
around the points marked with a cross show the systematic uncertainties in the Au-|-Au measurements. 



FIG. 15. Scaled third-order Fourier coefficients V 3 {pt)/( ssN^^^^) for charged hadrons measured at midrapidity in Gu-|-Au and 
Au-|-Au collisions at = 200 GeV. In each panel, the vzipr) values measured in the centrality classes marked in the 

figure, are scaled by the average third-order participant eccentricity £3 in the initial state of the collisions as determined by 
a MC Glauber calculation described in the text, and the corresponding number of nucleon participants A^p^rt- The symbols 

represent the scaled V 3 {pt)/[ ssN^^j.^) values, and the error bars show the statistical uncertainties. The shaded boxes indicate 
the systematic uncertainties in the Cu-|-Au measurements, and the fines around the points marked with a cross show the 
systematic uncertainties in the Au-|-Au measurements. 
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FIG. 16. The second-order Fourier coefficients U 2 (pt) for charged hadrons measured at midrapidity in Cu+Au collisions at 
= 200 GeV in comparison to hydrodynamics calculations for the centrality classes marked in each panel. The symbols 
represent the measured V 2 {pt) values, the error bars show the statistical uncertainties, and the shaded boxes indicate the 
systematic uncertainties. The theoretical calculations, shown with the solid and dashed lines, are performed with two different 
values of the specific viscosity rj/s marked in the figure. 
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FIG. 17. The third-order Fourier coefficients V 3 (pt} for charged hadrons measured at midrapidity in Gu-|-Au collisions at 
= 200 GeV in comparison to hydrodynamics calculations for the centrality classes marked in each panel. The symbols 
represent the measured V 3 (pt) values, the error bars show the statistical uncertainties, and the shaded boxes indicate the 
systematic uncertainties. The theoretical calculations, shown with the solid and dashed lines, are performed with two different 
values of the specific viscosity r)/s marked in the figure. 
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FIG. 18. The second-order Fourier coefficients U 2 (pt) for charged hadrons measured at midrapidity in Cu+Au collisions at 
= 200 GeV in comparison to AMPT model calculation for the centrality classes marked in each panel. The symbols 
represent the measured V 2 {pt) values, the error bars show the statistical uncertainties, and the shaded boxes indicate the 
systematic uncertainties. For the theoretical calculations, which are shown with fines, only statistical uncertainties are plotted. 
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FIG. 19. The third-order Fourier coefficients V 3 {pt) for charged hadrons measured at midrapidity in Gu-|-Au collisions at 
y'Sjyjy = 200 GeV in comparison to AMPT model calculation for the centrality classes marked in each panel. The symbols 
represent the measured V 3 (pt) values, the error bars show the statistical uncertainties, and the shaded boxes indicate the 
systematic uncertainties. For the theoretical calculations, which are shown with fines, only statistical uncertainties are plotted. 
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FIG. 20. The top panels show the maximum-normalized distribution of hrst-order participant plane angle Vlti computed 
from the initial coordinates of participant nucleons determined with (a) black disk, and (b) gray disk Monte Carlo Glauber 
simulations. The bottom panels show AMPT wi computed with respect to the impact parameter b, and ^<1 using (c) black disk 
and (d) gray disk Monte Carlo for the initial conditions. 
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We note that it is currently unknown whether the spec¬ 
tator neutrons bend toward or away from the interaction 
overlap region between the nuclei, and whether they are 
oriented along the impact parameter vector 6 , along the 
vector dti determined by the initial energy density in the 
overlap region, or some other vector. In fact, it is con¬ 
ceivable that spectators very close to the overlap region 
have a different behavior from spectators far away from 
the overlap. These ambiguities need resolution before a 
more direct theory to data comparison can be made. 

IV. SUMMARY 

Anisotropic flow coefficients for inclusive charged par¬ 
ticles and identified hadrons tt^ , , p, and p produced 

in Cu-I-Au collisions at = 200 GeV have been mea¬ 

sured by the PHENIX experiment at RHIC using event 
plane techniques. The vi, V 2 , and measurements were 
performed at midrapidity as a function of transverse mo¬ 
mentum pt over a broad range of collision centralities. 
Mass ordering was observed for low px in the identified 
particle measurements, as predicted by hydrodynamics. 

A system size comparison was performed for the in¬ 
clusive charged particles using previous PHENIX mea¬ 
surements at = 200 GeV of V 2 {pt) in Gu-|-Cu and 

Au-|-Au collisions, and V 3 {pt) in Au-|-Au collisions. The 
elliptic and triangular flow measurements between differ¬ 
ent systems and centrality selections were found to scale 
with the product of the initial participant eccentricity 
and the third root of the number of nucleon participants 
e„Vp'rf The system size dependence of the V 2 ,{pt) val¬ 
ues could also be described by participant eccentricity £3 
scaling alone. 

The inclusive charged-particle measurements were 
compared to theoretical predictions. In the vi measure¬ 
ment, we observed negative values at highp^, indicating 
that hadrons are emitted in the transverse plane pref¬ 
erentially in the hemisphere of the spectators from the 
Au nucleus, assuming that they moved outward from the 
interaction region and are aligned with the impact pa¬ 
rameter vector. The AMPT transport model calculations 
were found to be in agreement with the magnitude of the 
measured vi{pt) signals, but having the opposite sign. 
At low Pt (< 1 GeV/c) AMPT provides a reasonable de¬ 
scription of the triangular flow in all measured centrality 


classes that cover the 0%-30% range, and the elliptic flow 
measurements in the 0%-60% range. Event-by-event hy¬ 
drodynamics calculations with specific viscosity in the 
range rj/s = 0.08 — 0.16 reproduce the measured V 2 (pt) 
and v^(jpt) values. 
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